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Recent experiments with heavy ions and planned experiments with ultra intense lasers require 
nonperturbative solutions to quantum field theory for predicting and interpreting the results. To 
propel this theoretical direction, we solve the nonperturbative problem of an electron in a strong 
transverse confining potential using Hamiltonian light-front quantum field theory. We evaluate both 
the invariant mass spectra and the anomalous magnetic moment of the lowest state for this two-scale 
system. The weak external field limit of the anomalous magnetic moment agrees with the result of 
QED perturbation theory within the anticipated accuracy. 
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PACS numbers: ll.10.Ef, 11.15.Tk, 12.20.Ds 

Recent intense interest in strong-field dynamics, rang- 
ing from the anomalous enhancement of lepton produc- 
tion in ultrarelativistic collisions between heavy nuclei 
at RHIC [l[ and predicted photon yield depletion at the 
LHC [2|, to proposals for producing supercritical fields 
with next-generation laser facilities [3, |j| , points to the 
importance of developing new methods for solving QED 
in its nonperturbative domain. An ideal tool for such 
problems is Hamiltonian light-front formalism [see e.g. 
Ref.@|], in which the gauge theory is quantized on the 
light-front and the physical states are expanded in a 
Fock-space basis developed from the constituents. The 
Hamiltonian is represented as an operator acting on these 
Fock states. Since time is set along the light-front, the 
ground state of the free theory is also a ground state of 
the full interacting theory and the formalism is Lorentz 
frame independent. 

The light-front Hamiltonian approach provides a real- 
ization of Feynman's covariant "parton" model where the 
partons are the elementary fields used to define the Fock- 
space basis. It also has the appearance of a standard 
quantum many-body problem with the necessary quan- 
tum field theory features such as pair creation and an- 
nihilation. Diagonalization of the Hamiltonian provides 
amplitudes for evaluating experimental observables that 
are nonperturbative and relativistically invariant such as 
masses, form factors, and structure functions. The Fock- 
space dependence of observables is to be evaluated and 
one seeks to eliminate such dependence. 

We address the problem of an electron in a transverse 
harmonic cavity and solve for its mass spectra and other 
observables as a function of the external field strength 
over a range spanning the electron mass. To accom- 
plish this, we evaluate the QED Hamiltonian in light- 
front gauge on the light-front in a Fock-space consisting 
of electron states and electron plus photon states. We 
add the harmonic oscillator potential in the transverse 
direction to confine the system in those directions. We 
then solve for the eigenvalues, eigenvectors, and anoma- 
lous magnetic moment. The nonperturbative analysis 
presented in this Letter could be applicable to measure- 



ments of the (gyromagnetic) ratio of the spin precession 
to Larmor frequencies of a trapped electron in strong ex- 
ternal electromagnetic fields or intense time-dependent 
laser fields. This research also serves as a foundation for 
solving other quantum field theories at strong coupling, 
such as the light-front QCD Hamiltonian in the nonper- 
turbative domain. 

The question arises on how to implement a consistent 
renormalization program. Our specific choice is defined 
below. With our limits on the Fock-space and adopted 
renormalization scheme, we can already demonstrate the 
effects of mass renormalization but not yet coupling con- 
stant renormalization. A full renormalization program 
will ensue when we enlarge our Fock-space basis to in- 
clude electron-positron pairs. Such pairs produced at 
RHIC experience strong time-dependent EM fields. The 
total charge Z to t^i = Z± + Z2 of the two colliding nu- 
clei can exceed 137, indicating the need to treat the pro- 
duction and propagation properties with nonperturba- 
tive methods. In particular, the renormalization scale 
entering the running QED coupling is set by the pho- 
ton virtuality which would be expected to be of order 
■Ztotaio^e- The renormalization scale p 2 appearing in the 
electron running mass m(p 2 ) is of similar order. Such a 
full renormalization program will also allow comparison 
with higher order perturbative calculations such as those 
for an electron in a Penning trap [|| . 



We define our light-front coordinates as s J 



x u ±x J 



x 1 - = (x 1 ^ 2 ), where the variable x + is light-front time 
and x~ is the longitudinal coordinate. We adopt x + = 0, 
the "null plane" , for our quantization surface. Here we 
adopt basis states for each constituent that consist of 
transverse 2D harmonic oscillator (HO) states combined 
with discretized longitudinal modes, plane waves, satis- 
fying selected boundary conditions. This basis function 
approach follows Q and is supported by successful anti- 
de Sitter-QCD models 0. 

The HO states are characterized by a principal quan- 
tum number n, orbital quantum number m, and HO en- 
ergy f2. Working in momentum space, it is convenient to 
write the 2D oscillator as a function of the dimensionless 
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variable p = \p \J ^/ MqQ,, and Mq has units of mass. The 
orthonormalized HO wave functions in polar coordinates 
(p, ip) are then given in terms of the generalized Laguerre 
polynomials, L^(p 2 ), by 



$nm(P,<P) = (p<p\nm) 
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with eigenvalues E n>m = (2n + \m\ + The HO wave- 
functions have the same analytic structure in both coor- 
dinate and momentum space, a feature reminiscent of a 
plane- wave basis. 

The longitudinal modes, ipk, in our basis are defined 
for — L < x ~ < L with periodic boundary conditions for 
the photon and antiperiodic boundary conditions for the 
electron: 



ipk{x ) = 



1 



2L 



e L 



k x 



(2) 



where k — 1, 2, 3, ... for periodic boundary conditions (we 
neglect the zero mode) and k — |, |, |, ... for antiperiodic 
boundary conditions. The full 3D single-particle basis 
state is defined by the product form 

^k,n,m( X ~iPi p) = 1pk(x^)^ n , m (p,p)- (3) 

Following Ref. [I| we introduce the total invariant mass- 
squared M 2 for the low-lying physical states in terms of 
a Hamiltonian H times a dimensionless integer for the 
total light-front momentum K 



M 2 + Pi Pi -» M 2 + const = P + P~ 



KH 



(4) 



where we absorb the constant into M 2 . For simplicity, 
the transverse functions for both the electron and the 
photon are taken as eigenmodes of the trap. The non- 
interacting Hamiltonian Ho = 2MqP~ for this system is 
then defined by the sum of the occupied modes i in each 
many-parton state: 



Hn = 



2M a n 



K 



E 



2m + \rrii 



1 + m?/(2M n) 



,(5) 



where fhi is the mass of the parton i. The photon mass 
is set to zero throughout this work and the electron mass 
fh e is set at the physical mass 0.511 MeV in our non- 
renormalized calculations. We also set Mq = m e . 

The light-front QED Hamiltonian interaction terms we 
need are the electron to electron-photon vertex, given as 



, (6) 
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and the instantaneous electron-photon interaction, 



a 2 f — t+ 
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(7) 
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FIG. 1: (color online). Eigenvalues (multiplied by K) for a 
nonrenormalized light-front QED Hamiltonian which includes 
the electron-photon vertex and the instantaneous electron- 
photon interaction without counterterms. The cutoffs for the 
basis space dimensions are selected such that K increases si- 
multaneously with the Nmax- 



where the coupling constant g 2 = Aira, and a is the 



fine structure constant taken to be a 



in this 



137.036 

work. The nonspinflip vertex terms of Eq.(|6]) are oc MqQ, 
whereas spinflip terms are oc \fMo£lm e . Selecting the ini- 
tial state electron helicity in the single electron sector al- 
ways as "up" the process e — > ej is nonzero for three out 
of eight helicity combinations, and the process ej — > 67 
is nonzero only with all four spin projections aligned (two 
out of 16 combinations), resulting in a sparse matrix. 

We implement a symmetry constraint for the basis 
by fixing the total angular momentum projection J z = 
M + S = |, where M = J2 t 

rrii is the total azimuthal 
quantum number, and S — J2i s i * ne total spin projec- 
tion along the x~ direction. For cutoffs, we select the 
total light-front momentum, K, and the maximum to- 
tal quanta allowed in the transverse mode of each one or 
two-parton state, N max , such that 



i i 

^2 2n i + \ m i \ + 1 < N max , 



(8) 
(9) 



where, for example, ki defines the longitudinal modes of 
Eq.© for the i th parton. Equation ((HJ signifies total 
light-front momentum conservation written in terms of 
boost-invariant momentum fractions, a;,-. Since we em- 
ploy a mix of boundary conditions and all states have 
half-integer total K, we will quote K values rounded 
downwards for convenience, except when the precise 
value is required. 

In Fig[T]we show the eigenvalues (multiplied by K) for 
a nonrenormalized light-front QED Hamiltonian given 
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in Egs. ([gl6l7]) . with fixed Q = 0.05 MeV and simulta- 
neously increasing K and N max . The resulting dimen- 
sion of the Hamiltonian matrix increases rapidly. For 
N max = K = 2, 10, and 20, the dimensions of the corre- 
sponding symmetric d x d matrices are d = 2, 1670, and 
26 990, respectively. 

The number of the single electron basis states, consid- 
ering all the symmetries, increases slowly with increasing 
N max = K cutoff. For N max = K = 2, 10, and 20 the 
number of single electron basis states is 1, 5, and 10, re- 
spectively. Our lowest-lying eigenvalue corresponds to a 
solution dominated by the electron with n — m = 0. The 
ordering of excited states, due to significant interaction 
mixing, does not always follow the highly degenerate un- 
perturbed spectrum of Eq. ((5|) . States dominated by spin- 
flipped electron-photon components are evident in the 
solutions. Nevertheless, the lowest-lying eigenvalues ap- 
pear with nearly harmonic separations in Figfl] as would 
be expected at the coupling of QED. The multiplicity 
of the higher eigenstates increases rapidly with increas- 
ing N max = K and the states exhibit stronger mixing 
with other states than the lowest-lying states. In prin- 
ciple, the electron-photon basis states interact directly 
with each other in leading order through the instanta- 
neous electron-photon interaction, but numerically the 
effect of this interaction is very weak, and thus does not 
contribute significantly to the mixing. Even though we 
work within a Fock-space approach, our numerical results 
should approximate the lowest order perturbative QED 
results for sufficiently weak external field. 

In Figj2] we show the results for the square root of the 
electron anomalous magnetic moment (scaled), J S[i/g 2 , 
as a function of £2 obtained from the lowest mass eigen- 
state. That is, we plot the magnitude of the proba- 
bility amplitude that electron has its spin flipped rela- 
tive to the single electron Fock-space component in the 
range where the results are converged. Since our system 
is in an external field, the lowest physical mass eigen- 
state (not known experimentally) can deviate from the 
free electron mass. Therefore, without renormalization, 
we only consider cases where the mass eigenvalue falls 
within 25% of the free electron mass. At zero external 
field we may compare our 5 fj, to the QED one- loop con- 
tribution to the electron anomalous magnetic moment, 
the Schwinger result £- [9J. That is we compare our re- 
sults with ^a/ {2itg 2 ) = ^/1/8tt 2 . For even N max = K 
the results converge rapidly for N max = K > 14. The 
results for odd cutoffs (not shown) track even cutoff re- 
sults as N max — K increases. Below £)<0.05 MeV all 
interactions are quenched at fixed N max — K, and not 
converged, due in part to our requirement that the HO 
basis tracks the external field. 

Figure [2] also shows an extrapolation of the above re- 
sults for N max = K = 12, . . . , 20 to the zero exter- 
nal field limit £7 — Mev. We have only included 
the points above the peak at £7 > 0.05 MeV, where we 
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FIG. 2: (color online). Square root of the (scaled) electron 
anomalous magnetic moment as a function of the transverse 
external field for a sequence of increasing basis spaces (solid 
grey lines). These are nonrenormalized results where the mass 
eigenvalue falls within 25% of the free electron mass. The 
theoretical one-loop QED prediction ("Schwinger") result of 
0.1125, appropriate to SI = MeV, is indicated. The black 
solid line is a fit to the results for N„ lax = K = 12, . . . 20, 
points included into the fit are indicated by the markers in 
the legend. Extrapolation to zero external field yields 0.1121. 



have reasonable convergence. An excellent agreement 
with the results is obtained by a fit function /(O) = 
a(l+bfl 2 +cn 4 ) exp(-cZfi), with a = 0.1121. This is < 1% 
deviation from the Schwinger result of 0.1125, which is 
reasonable in light of our numerical accuracy and extrap- 
olation uncertainties. If we perform individual extrapo- 
lations for all the N max = K = 12, . . . , 20 results with 

0. 1 < fl < 1.4 MeV, a range spanning the electron mass 
scale, we obtain excellent fits with 0.1109 < a < 0.1134, 

1. e., remaining within 1.5% of the Schwinger result. 

In Ref.[3| we discussed possible divergences present in 
our framework, and anticipated a straightforward man- 
agement of the identified divergences. Here we renor- 
malize our results by applying a sector-dependent nor- 
malization scheme from Ref.[10|. In our present lim- 
ited Fock-space, we need only the mass counterterm Sm e . 
This 5m e is added to the mass term in the diagonal one- 
electron part of the Hamiltonian Eq.©. In the absence 
of a known experimental mass for renormalization due 
to the external field, we adjust 5m e such that the low- 
est eigenstate remains at KEq = m 2 + Mrjfi. That is, 
we simply adopt the free electron mass for the renormal- 
ized mass, and keep the coupling constant g 2 unchanged. 
We emphasize that our choice for the renormalized mass 
and for the coupling constant are, in principle, valid for 
the case of zero external field only. Measurements for 
electrons in a trap (see, e.g., 



11]) could provide results 



leading to more precise renormalization parameters, but 
this aspect is beyond the scope of this Letter. 
Fig J3] we present ^S^i/g 2 for N max 
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FIG. 3: (color online). Individual fits (solid black lines) to the 
renormalized results for square root of the (scaled) electron 
anomalous magnetic moment for N max = K = 10, ... , 20. 
The inset shows the continuum limit extrapolation of the zero 
external field results in the main panel as a function of 1/K. 



10, . . . , 20 from the renormalized QED Hamiltonian of 
Eq.©, with Sm e , and Eqs. (|6l7[) . To eliminate possi- 
ble effects from the peak at Q ~ 0.05 MeV in FigH 
we only include results with the external field S7 > 0.2 
MeV. Again, individual fits of the form /(fi) = a(l + 
bfl 2 + cSl A ) exp(— dSl) are an excellent representation of 
our results. The range of the extrapolated values is 
0.1077 < a < 0.1216. 

The convergence with an increasing cutoff is now less 
rapid than in the nonrenormalized case shown in Figl^l 
In order to approach the continuum limit N max — K — > 
oo, we perform further extrapolation to the zero-57 re- 
sults of Fig|3l The inset of Figj3] shows linear extrap- 
olation of the results of the main figure in 1/K to the 
continuum limit N max — K — > oo. To verify the stabil- 
ity of the results, an extrapolation based on the f2 > 0.1 
MeV fits (not shown) is also given. The extrapolated 
continuum values are 0.1362 (0.1383) for VL > 0.2(0.1), 
respectively, and thus about 20% above the Schwinger re- 
sult 0.1125. An enhancement of this magnitude was also 
observed in related works, Refs. 



12J,113| and Refs.[LJ,[l 



where the one-photon truncated light-front Hamiltonian 
was regulated with Pauli-Villars regularization scheme. 
With Pauli-Villars regularization as well as in our renor- 
malized results, interpreted from a perturbation theory 
perspective, the intermediate state propagators are devel- 
oped from a dynamical (nonperturbative) electron mass 
rather than using the unperturbed mass needed for direct 
comparison with perturbation theory. 

In our approach, the HO parameters £7, Mq, the elec- 
tron mass m e , and the total longitudinal momentum 
K appear as prefactors for the matrix elements in the 
Hamiltonian. Therefore, we can rather straightforwardly 
vary the size of the Hamiltonian matrix by keeping N max 



fixed, and changing K alone. We studied the continuum 
limit of y/ Sfi/g 2 by setting N max = 20 and increasing 
K in units of 10, from K = 10 to K = 50. The di- 
mension of the Hamiltonian matrix then increases from 
d = 11790 to d = 69590. The extrapolated results range 
between 0.1148 < a < 0.1259, and show a good con- 
vergence pattern. Linear extrapolation of these results, 
analogously to Fig(3j to the continuum limit K — > oo are 
0.1288 (0.1290) with il > 0.2(0.1) MeV, - 15% above 
the Schwinger value. 

In summary, we have evaluated properties of an elec- 
tron in a nonperturbative external harmonic oscillator 
potential. We have taken the weak external field limit of 
the electron anomalous magnetic moment, and obtained 
results compatible with QED perturbation theory with 
reasonable accuracy. Our framework can be extended 
by incorporating higher Fock-space sectors and adopting 
external strong fields relevant to heavy ion collisions and 
to future high-intensity laser facilities. Applications to 
QCD will proceed with the adoption of recently devel- 
oped color-singlet basis enumeration techniques [7[. 
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